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Degenerate sterile neutrino dark matter in the cores of galaxies 
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Abstract. We study the distribution of fermionic dark matter at the center of galaxies using NFW, Moore and isothermal 
density profiles and show that dark matter becomes degenerate for particle masses of a few keV and for distances less than 
a few parsec from the center of our galaxy. A compact degenerate core forms after galaxy merging and boosts the growth of 
supermassive black holes at the center of galaxies. To explain the galactic center black hole of mass of ~ 3.5 x 10 6 M o and a 
supermassive black hole of ~ 3 x 10' M Q at a redshift of 6.41 in SDSS quasars, we require a degenerate core of mass between 
3 x 10 3 M o and 3.5 x 10 6 M o . This constrains the mass of the dark matter particle between 0.6 keV and 82 keV. The lower limit 
on the dark matter mass is improved to 7 keV if exact solutions of Poisson's equation are used in the isothermal power law 
case. We argue that the constrained particle could be the long sought dark matter of the Universe that is interpreted here as a 
sterile neutrino. 
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1. Introduction 

Precision observations of the cosmic microwave background 
and of large scale structure confirm the picture in which 96% 
of the matter density of the Universe is made of dark energy 
and dark matter (DM), which could be revealed only by their 
gravitational interaction (Spergel et al. 2003 1. The nature of 
these two forms of matter is still unknown. Many candidates 
have been proposed for DM. These include Cold Dark Matter 
(CDM) particles of masses heavier than lGeV/c 2 (Bertone et 
al. I2005i and Warm Dark Matter (WDM) particles such as 
sterile neutrinos (Dodelson & Widrow ri994> . Recently, there 
has been a renewed interest in sterile neutrinos as candidates 
for DM as they could be a natural extension of the minimal 
standard model (MSM) of electroweak interactions (Weinberg 
I1967I Glashow ll961> . One of these sterile neutrinos could be 
the DM, while the other two could help explain baryogenesis 
(Asaka & Shaposhnikov 2005 Shaposhnikov|2UU(>J. 

Numerical cosmological N-body simulations suggest DM 
density profiles which follow p ~ r r law, with y « 3 in the 

outer parts of the halos and 1 ~ y ~ 2 inside a few kpc (Navarro 
et al. fl997l Moor e et al. IT9981 IT9991 Klypin et al. 2002] and 
Power et al. 2003 1. Although the Milky Way is well studied in 
the inner 3-10 kpc region, little is known about the DM distri- 
bution on smaller scales, i.e. r ~ 0.1 pc where the gravitational 



potential of baryons dominates over DM. In addition, there is a 
mounting evidence that a black hole of mass M ~ 3.5 1O 6 M 
dominates the mass distribution in the inner one parsec of the 
Galaxy (Schodel et al. l2003l Ghez et al. 12005V Thus the inves- 
tigation of the distribution of DM around the black hole in the 
Galactic center is important as it could give some hints on the 
nature of the DM particle (Gnedin & Primack 2005 Bertone & 
Merritt 2005 ). Recently, we have established that the inner DM 
density profile scales as p ~ r~ 3 ^ 2 under the assumption that the 
black hole feeds from degenerate fermionic DM (Munyaneza 
& Biermann 2005 1. It was then shown that a DM particle of or- 
der of m s ~ 10 keV could explain the growth of supermassive 
black holes of 10 6 to 1O 9 M from stellar seed black holes. It 
is interesting to note that a sterile neutrino mass in overlapping 
mass range, with a small mixing angle with the active neutri- 
nos, has also been suggested by Kusenko (20041 to explain the 
high velocities of pulsars at birth in supernovae. Moreover, a 
DM particle mass in the keV range can be constrained from 
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X-ray background studies (Drees & Wright 2000, Abazajian et 
al. HSU and Dolgov & Hansen 2002). Here, it is worth men- 
tioning that an upper limit of a sterile neutrino mass of 6.3 keV 
from X-ray background seems to be in conflict with the lower 
limit of 14 keV from SDSS Lyman-alpha forest (Abazajian & 
Koushiappas 2006 1. However, the non-thermal phase distribu- 
tion of the DM particles has the potential to modify the Lyman 
alpha limits, while leaving the X-ray limits intact, so elim- 
inating the contradiction with this additional degree of free- 
dom. Depending on the specific model for the production of the 
DM particles, their initial phase space distribution is possibly 
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largely sub-thermal. Recently, Biermann & Kusenko ( 2006|l es- 
tablished that the decay of a such sterile neutrino could speed 
up the formation of molecular hydrogen and boost the early 
star formation and reionisation in agreement with the WMAP 
3-year results (Spergel et al. 2006 1. 

The purpose of this Letter is to study the constraints on the 
DM particle mass in the central region of the galaxy. Given that 
there is still uncertainty about the DM profile around the cen- 
tral black hole, we will assume standard Navarro-Frenk- White 
(NFW), Moore or isothermal gas sphere profiles and investi- 
gate what happens if the DM particles become degenerate and 
form a degenerate dark matter star also called a fermion ball 
(Munyaneza & Viollier 2002 Munyaneza & Biermann I2005> . 
We then get the lower limits on the DM particle mass by assum- 
ing that the mass of the fermion ball cannot be more than the 
mass of 3.5 x 1O 6 M for the Galactic center black hole. An up- 
per limit to the DM mass will be established from the condition 
that the mass of the fermion ball should be greater than about 
3 x 10 3 M o . The choice of this mass comes from the growth 
mechanism arguments to form supermassive black holes of 
mass of 10 9 M Q in SDSS quasars (Munyaneza & Biermann 
2005 1. Moreover, a mass of 3 x 1O 3 M is the upper limit of 
masses for which there is no black hole in galaxies such as M33 
(Gebhardt et al. l2"0Tm & Barth et al.|2005). Thus the Galactic 
center provides us with a fertile testing ground for DM theories 
and cosmological structure evolution. 

2. Degenerate cores 

We study the following density profiles obtained through 
galaxy mergers 



P(r) = po 



1 



1 



(1) 



with s < 3 and p a is the characteristic density at a size of ro = 
3 kpc. The value of po will be found from the constraint that the 
DM mass enclosed at r — 3 kpc is about 1O 1O M for the Milky 
Way galaxy (Klypin et al. 2002i. The index s characterizes the 
density power law slope and we will use the value of s = 1 
for NFW pro file (Navarro , Frenk & White IT9971 . s = 1.5 for 
Moore et al. ( 119981119991 profile and s = 2 for the isothermal 
gas sphere case. The full distribution function corresponding to 
the the density of equation Q can be found in Hansen et al. 
(2005 1. Of course, this distribution function does not allow for 
the Fermi Dirac nature of the DM particles. The mass enclosed 

within a radius r and for small r ~ ro = 3 kpc is given by 



M(r) 



= I Anr p(r)dr - r, 
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The rotational velocity v ror corresponding to the mass enclosed 
M(r) is 
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Here, we will consider only DM contribution to v rot . The nu- 
merical simulation results of inner density power law density 



distributions suggest that for a large number of particles, their 
kinetic energy is strongly reduced and the more so, the deeper 
they are near the center. The central particle populations are 
essentially cooled in the core in the violent relaxation of the 
merger, with a local temperature T 
condition can be written as 



r s . The degeneracy 
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Here, m s is the DM particle mass to be constrained, n = p/m s is 
the DM particle number density, ti is the Planck's constant and 
gf is the spin degree of freedom, i.e. gf = 2 for Majorana and 
gf = 4 for Dirac's particles. The constraint 10} implies a degen- 
erate Fermi Dirac distribution which of course does not have a 
long power law high energy tail. The Fermi Dirac distribution 
function /fd is given by 



f FD (E) = K- 



1 
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where E is the DM total energy, p. the chemical potential, T the 
temperature and K is the normalisation constant and N the total 
number of DM particles. For high energies, i.e. E - p » kT, 
the Fermi - Dirac distribution function becomes Maxwellian, 
i.e. foE ~ e~ (E ~v)/i<T w hereas direct numerical simulations 
favour a velocity distribution with a long power law high en- 
ergy tail (Hansen & Stadel l2006l Hansen et al. 2005 1 with the 
following profile 



/(v) = 



l-(l-q)\ 

1 K\(T 



(6) 



where q and k\ are free parameters. The latter distribution be- 
comes Maxwellian for q — 1 . We also note that for v — > 0, the 
distribution function given by equation (|6} becomes indistin- 
guishable from a degenerate Fermi Dirac distribution at abso- 
lute zero Temperature. To reproduce the characteristic features 
of the distribution given by equation (|5Jl in the Fermi Dirac dis- 
tribution, we introduce a multi-temperature distribution func- 
tion as follows 



f(E)=KY 



1 



^ exp(^) + 1 



(7) 



where T„ = 2" To and the chemical potential p„ is choosen from 
the condition that exp {j^rf^) — pr- In result we get 



p„ = -kTon2" In p, 



(8) 



where n is the power index and p" shows by how much the 
distribution function has been lowered. In the limit of n — > oo, 
the sum from equation can be converted into the following 
integral 



f(E) 



1 



dT 



to exp(^) + 1 T 



(9) 



Thus, the chemical potential and the temperature parameters 
allow us to achieve a high energy power law tail in the Fermi 
Dirac distribution. For realistic velocity distributions obtained 
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from large scale simulations such as given by equation (j5Jl, the 
dominant contribution to the number of particles is still the ba- 
sic temperature given by /qcr, demonstrating that a Fermi Dirac 
distribution is a viable first approximation. The modification to 
the degeneracy criterion of equation (0} is small. In addition, 
the density powerlaw slope s in NFW, Moore and isothermal 
distributions can be related to our parameter p so that the num- 
ber of particles at the temperature Tq dominates by far over all 
other partial sums. 

The characteristic Fermi velocity is denoted by v and will 
be taken here as the rotational velocity, v ro ,. Using equations 
and 0, we find that the degeneracy condition @ becomes 
fulfilled for particle masses of a few keV and for distances less 
than a few parsec. Here, we point out that the physics of de- 
generate fermion balls that obey equation ® was studied in a 
series of papers (Bilic et al. 1999 Tsiklauri & Viollier 1998 
Munyaneza et al. [T998|[T999l Munyaneza & Viollier l2002t . 

Introducing the velocity dispersion <x 2 = -=-&, where Mo = 

T° Po is the total mass of DM contained within a size of 

i-S 

ro = 3 kpc. As both the velocity v and the number density n 
depend on the radial coordinate, the degeneracy condition 10} 
allows to get a constraint on the DM particle mass m s and the 
radius r. Let us denote by R dm the radius at which DM be- 
comes degenerate 



Rt 



I 8f 



4n 



Using the last expression for the radius, we can calculate the 
mass M DM of DM enclosed within the above radius R DM to be 



M nM = 



\ 8f 



An 



V c I \ ro / \m s 



where m v \ = (fic/G) 1 ^ 2 and l p \ = (/zG/c 3 ) are the Planck's 
mass and length, respectively. We wish to emphasize that this 
is just the mass derived from the original power law density 
distribution, and so does not include any additional DM accre- 
tion to the DM degenerate configuration. We then rewrite the 
equations (I10> . il It for each value of the steepness index s. 
Thus, for NFW profile, i.e. s = 1 we get 



M DM = 5.24 m, 



a\ 6 ' 5 //pA 3/5 /mpi 
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(12) 



(13) 



for the total mass and the size of the degenerate core, respec- 
tively. For Moore profile, i.e. s = 1.5, the total mass Mdm and 
the size R dm of the degenerate core are given by 



(14) 



(15) 




[keV] 



Fig. 1. The total mass Mdm of the fermion ball as a function 
of the fermion mass m s . The total mass Mdm scales with the 
the DM particle mass m s as m^ 16 ^ 5 , m, and as m~ 2 for NFW, 
Moore and isothermal power law, respectively. Two horizontal 
lines at M D m = 3.5 x 1O 6 M and M DM = 3 x 10 3 M Q have been 
drawn to get the lower and upper limits on the mass of the DM 
particle. 



and finally for the isothermal gas sphere, i.e for s = 2, the mass 
Mdm and the size Rdm scale as 



M D m = 2.58m 



re' 



= 1.29/ 



■(16) 



In order to investigate the behaviour of the degenerate DM 
mass Mdm and its size Rdm as functions of the DM mass 
m s , we fix the total DM mass Mo = 1O 1O M at a distance of 
ro = 3 kpc and this gives a velocity dispersion of cr — 84.65 
km/s (Ferrarese & Merritt 2000 1 for the three density profiles. 
Here, only DM contribution to the total mass is taken into ac- 
count. 

' (l^We then assume that the total mass of the degenerate core 
should be between 3 x 10 3 M o and 3.5 x 1O 6 M . There are 
two arguments for the choice of this lower limit on the mass 
of fermion ball. First, this mass is needed to boost the growth 
of stellar seed black holes to ~ 1O 6 M which then grow to 
~ 1O 9 M supermassive black hole in SDSS quasars by bary- 
onic accretion (Munyaneza & Biermann 2005). Secondly, this 
is the threshold mass, below which no black hole will be 
formed as in the case of M33 (Gebhardt et al. 120011 & Barth et 
al. l2005Y 

For NFW profile, the mass of the DM particle should be in 

the range of 0.6 keV ~ m s ~ 6 keV. A range of 1 keV ~ 

m s ~ 14 keV is found for Moore profile and finally we obtain 

2 keV ~ m s ~ 82 keV for the isothermal gas sphere. These 
limits on the DM particle mass are illustrated in Fig. 1 where we 
plot the total mass of the degenerate core as a function of the 
DM mass m s . Using these lower limits for the DM particle, we 
predict a density of (1 - 100) M pc~ 3 in the case of NFW and 
Moore profiles at the size of a few parsec, which corresponds 
to the values obtained in Klypin et al. (2002). In the case of the 
isothermal gas sphere, we predict a much higher density of a 
few 10 5 M Q pc~ 3 at a size of about 1 pc. 

Here, we want to point out that the total mass of DM core 
in the isothermal gas sphere case is increased by a factor of 
nine if we use the full numerical solutions of Poisson's equa- 
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tion for a self gravitating system of fermions (Munyaneza & 
Biermann 12005b . Thus, in order to get a degenerate core of 
3.5 X 10 6 M o with a velocity dispersion of cr — 84.65 km/s, 
we get a fermion mass of m s * 7 keV instead of lower limit 
of m s ~ 2 keV. Using a fully developed distribution function 
such as the one given in equation would result in further 
changes of the obtained limits. However, we do not expect any 
major modifications since the number of particles in the tail is 
much smaller compared to the main body of the phase space 
distribution. 

3. Conclusion and Discussion 

We have determined the constraints on the DM mass under the 
assumption that the NFW, Moore or isothermal density profiles 
become degenerate near the center. Assuming that the mass of 
the degenerate fermion ball is between about 3 x 1O 3 M and 
3.5 x 10 6 M o , we have found that the mass of the DM parti- 
cle should have a lower limit of about 7 keV in the isothermal 
power law case. The above limit on the DM mass is obtained 
from a simple model which does not take into account baryonic 
matter and a proper model for DM with anisotropic distribution 
in phase space (Evans & An 2006 Hansen & Stadel 2006 ) and 
baryonic matter would modify these limits. The obtained mass 
range is in full agreement with the Tremaine-Gunn lower bound 
of about 0.5 keV on the mass of any fermionic DM when ap- 
plied to dwarf spheroidal galaxies (Tremaine & GunnQ]H3). 
The constraints obtained on DM particle mass overlap with the 
upper limit on sterile neutrino masses of about 8 keV to obey 
the X-ray emission constraints from the Virgo cluster observa- 
tions (Abazajian 2005). In addition, a sterile neutrino mass in 
the mass range of 2-20 keVs was derived to explain the high 
velocities up to 1000 km/s experienced by pulsars at birth in 
supernovae (Kusenko 2004 Barkovitch et al. 2004 and Fuller 
et al. 2003 ). Moreover, the obtained constraints agree well with 
the lower bound of m s ss 2 keV on the mass of sterile neutrinos 
obtained from the analysis of the Lyman- a forest data (Hansen 
et al. 127701 Viel et al. l2005l and Abazaiian 12005V Recently, 
Biermann & Kusenko ( 2006 1 have shown that the decay of such 
a sterile neutrino could help initiate star formation in the early 
Universe and the detection of such an X-ray line would con- 
firm whether such a sterile neutrino exists or not and this could 
be done by observing the X-ray decay line using XMM and 
Chandra satellites (Boyarski et al. 2006 Riemer-S0rensen et 
al. l2006Tl . 

To summarize, we find that NFW, Moore and isothermal 
density profiles become degenerate for fermion DM particle of 
keV masses at a size of the order parsec. Assuming that su- 
permassive black holes grow from degenerate fermion cores of 
masses of a few 1O 3 M to ~ 1O 6 M at galactic centers, we have 
found that the sterile neutrino mass should be in the range from 
0.6 keV to about 82 keV, with an improved lower limit of about 
7 keV in the isothermal power law density profile. As this range 
of sterile neutrino masses overlaps with other results on sterile 
neutrino masses discussed in this Letter, this leaves the sterile 
neutrino to be an excellent candidate for DM. The formation 
and growth details of degenerate fermion balls is beyond the 
scope of this Letter and will be discussed in another paper. 
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